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ABSTRACT 

The inner region of the accretion disk onto a rotating magnetized central star 
(neutron star, white dwarf or T Tauri star) is subjected to magnetic torques which 
induce warping and precession of the disk. The origin of these torques lies in the 
interaction between the surface current on the disk and the horizontal magnetic field 
(parallel to the disk) produced by the inclined magnetic dipole: The warping torque 
relies on the radial surface current generated by the twisting of the vertical field 
threading the disk, while the precessional torque relies on the azimuthal screening 
current due to the diamagnetic response of the disk. Under quite general conditions, 
there exists a magnetic warping instability in which the magnetic torque drives the 
disk plane away from the equatorial plane of the star toward a state where the disk 
normal vector is perpendicular to the spin axis. Viscous stress tends to suppress the 
warping instability at large radii, but the magnetic torque always dominates as the disk 
approaches the magnetosphere boundary. The magnetic torque also drives the tilted 
inner disk into retrograde precession (opposite to the rotation of the disk) around the 
stellar spin axis. Moreover, resonant magnetic forcing on the disk can occur which may 
affect the dynamics of the disk. 

The magnetically driven warping instability and precession may be related to a 
number of observational puzzles. Examples include: (1) Spin evolution of accreting 
X-ray pulsars: It is suggested that the observed torque reversal of the disk-fed 
magnetized neutron stars is associated with the wandering of the inner disk around 
the preferred perpendicular state. (2) Quasi-periodic oscillations in low-mass X-ray 
binaries: The magnetic torque induces disk tilt, making it possible to explain the 
observed low-frequency QPOs using disk precession. (3) Super-orbital periods in a 
number of X-ray binaries as a result of warped, precessing disks. (4) Photometric 
period variations of T Tauri stars. 



Subject headings: accretion disks - instabilities - stars: neutron - X-rays: stars - stars: 
T Tauri - stars: magnetic fields 
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1. 



INTRODUCTION 



Interaction between an accretion disk and a magnetized central object lies at the heart 
of the physics of a variety of astrophysical systems, including accreting neutron stars, white 
dwarfs and pre-main-sequence stars (e.g., Frank et al. 1992; Hartmann 1998). The basic picture 
of disk-magnetosphere interaction was first outlined by Pringle & Rees (1972) following the 
discovery of accretion-powered X-ray pulsars. These are rotating, highly magnetized {B ~ 10^^ G) 
neutron stars that accrete material from a companion star, either directly from a stellar wind, 
or in the form of an accretion disk. The strong magnetic field disrupts the accretion flow at 
the magnetospheric boundary (typically at a few hundreds neutron star radii), and channels the 
plasma onto the polar caps of the neutron star. The magnetosphere boundary is located where 
the magnetic and plasma stresses balance, 



rate, ry is a dimensionless constant of order unity (For an aligned dipole, estimate of rj ranges from 
0.5 to 1; e.g., Pringle & Rees 1972; Lamb, Pethick & Pines 1973; Ghosh & Lamb 1979a,b, 1992; 
Aly 1980; Arons 1993; Wang 1995). In low-mass X-ray binaries containing weakly magnetized 
(B ~ 10® G) neutron stars, the magnetosphere lies close to the stellar surface; the neutron stars 
are spun up and eventually become millisecond pulsars (Alpar et al. 1982; see Phinney & Kulkarni 
1994 or Bhattacharya 1995 for a review). Similar magnetosphere-disk interaction also occurs in T 
Tauri stars (e.g., Hartmann 1998), where the stellar magnetic field {B ~ 10^ G at the surface) 
strongly affects the accretion flow, and provides a key ingredient to explain the anomalous slow 
rotation of these objects (e.g., Konigl 1991; Cameron & Campbell 1993; Shu et al. 1994). Finally, 
models of magnetized accretion disks have been applied to explain the observed UV and X-ray 
emission from some intermediate polars (the DQ Her subclass of CVs) (e.g., Yi & Kenyon 1997). 

A large number of theoretical papers have been written on the subject of the interaction 
between accretion disk and a magnetized star (e.g., Pringle &: Rees 1972; Ghosh & Lamb 1979a,b; 
Aly 1980; Lipunov & Shakura 1980; Anzer & Borner 1980,1983; Arons 1987,1993; Wang 1987,1995; 
Aly & Kuijpers 1990; Spruit & Taam 1990,1993; Shu et al. 1994; van Ballegooijen 1994; Lovelace 
et al. 1995,1998; Li, Wickramasinge & Riidiger 1996; Campbell 1997), and numerical study of 
this problem is still in its infancy (e.g.. Stone & Norman 1994; Hayashi et al. 1996; Goodson 
et al. 1997; Miller & Stone 1997; see also Toropin et al. 1999). Beyond the simple notion of 
magnetosphere (with the Alfven radius estimated by eq. |1-1| ), there is little consensus on the 
range and strength of magnetic interaction outside the magnetosphere boundary, on the efficiency 
of magnetic field dissipation in/outside the disk, and on how or where the plasma attaches 
to the field lines. Outstanding issues include whether the disk excludes the stellar magnetic 
field by diamagnetic currents or the field can penetrate and thread a large fraction of the disk, 
whether the threaded field remains closed (connecting the star and the disk) or becomes open by 




(1-1) 



where M and n are the mass and magnetic moment of the central object, M is the mass accretion 
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differential shearing, and whether /how magnetically driven wind is launched from the disk or the 
magnetosphere/corotation boundary. 

This paper introduces some new physical effects associated with disk accretion onto rotating 
magnetized stars. Most previous papers have, for simplicity, assumed that the stellar magnetic 
dipole is aligned with the spin axis. Essentially all previous authors, with the exception of 
Lipunov & Shakura (1980) (see §2.1), have assumed that the spin is aligned with the disk angular 
momentum — There is a good reason for this, as it is natural to think that the star attains 
its spin angular momentum from the accretion disk. However, when these assumptions are 
abandoned, new physical effects are revealed: We show that under quite general conditions, the 
stellar magnetic field can induce disk warping and make the disk precess around the spin axis. As 
the accreting plasma approaches the magnetosphere, there is a tendency for the disk plane to be 
driven away from the equatorial plane of the star toward being aligned with the spin axis (§2 and 
§4) . The origin of the precessional torque and magnetic warping instability lies in the interaction 
between the surface current on the disk and the horizontal magnetic field (parallel to the disk) 
produced by the stellar dipole. The electric current on the disk is an inevitable consequence of 
magnetic-field - disk interaction, although its actual form is uncertain and depends on whether 
the disk is diamagnctic or has a vertical field threading it, thus on the dissipative processes in 
the disk and the magnetosphere (§2). The magnetic force on the disk also give rise to vertical 
and epicyclic resonances which may drive disk perturbations (§3). We show that the magnetically 
driven warping and precession of accretion disks can potentially explain a number of outstanding 
puzzles related to disk accretion onto magnetized stars (§6). 

This paper is organized as follows. In §2 we calculate the precessional and warping torques on 
the disk (or a ring of disk) due to the stellar magnetic field for several models based on different 
assumptions about magnetic field - disk interaction. The torques are different for different models, 
but they are always present and are of the same order of magnitude. The key physics of the 
torques are summarized at the beginning of §2. Resonances due to magnetic forces on the disk are 
discussed in §3. In §4 we study the dynamics of a warped disk under the influence of the magnetic 
torques derived in §2. The criterion for the magnetically driven warping instability (including the 
effect of viscosity) is derived, and a magnetic "Bardeen-Petterson" effect is discussed. Section 5 
addresses the question of how the spin of the central star is affected by the warped, precessing 
disk. In §6 we discuss/speculate several astrophysical applications of our theory, including the 
spin evolution of accretion-powered X-ray pulsars, quasi-periodic oscillations in low-mass X-ray 
binaries, long-term periodic cycles in several X-ray binaries, and the variability and rotation of T 
Tauri stars. We conclude in §7 by discussing possible future studies along the line initiated in this 
paper. 

Because of the intrinsic uncertainties associated with the nature of magnetic field - disk 
interaction, in the main text we shall focus on generic features and rely on parametrized models. 
In Appendix A we discuss a global magnetized disk model. In Appendix B wc consider the issue 
of calculating the magnetosphere radius for nontrivial magnetic field and spin geometry. The 
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dynamics of an extreme form of disk (consisting of diamagnetic blobs) is examined in Appendix C. 

2. MAGNETIC PRECESSIONAL AND WARPING TORQUES 

In this section we calculate the magnetic torque on the disk (or a ring of disk) surrounding a 

rotating magnetic dipole. Because of the uncertainties on how the magnetic field behaves in the 
presence of an accretion disk (see §1), we shall consider several different models, some representing 
extreme situations. Our most "realistic" model is given in §2.3. 

The basic set-up used in our calculation is shown in Fig. 1. The disk angular momentum axis 
I is inclined by an angle /? with respect to the spin axis a). The stellar dipole momentum jl rotates 
around the spin axis oj. The angle of obliquity between the magnetic moment and the rotation 
axis is 6. 

The key physics responsible for the magnetic torques on the disk is as follows: 

(i) If the disk is diamagnetic (perfectly conducting) so that the vertical stellar field 
(perpendicular to the disk) cannot penetrate, an azimuthal screening surface current will be 
induced in the disk. This current interacts with the radial magnetic field from the star, and a 
vertical magnetic force results. While the mean force (averaging over the azimuthal) is zero, the 
uneven distribution of the force induces a net torque acting on the ring, making it precess around 
the spin axis. (For a nonrotating star, the disk will precess around the magnetic axis.) (§2.1). 

(ii) If there is a vertical field threading the disk (this either comes from the star or is 
carried intrinsically by the disk and detached from the star) , it will be twisted by the disk rotation 
to produce a discontinuous azimuthal field AB^ across the disk surface, and a radial surface 
current results. The interaction between this current and the stellar i?^ (which is not affected 
by the disk) gives rise to a vertical force, and the resulting torque tends to misalign the angular 
momentum of the ring with the spin axis (although in some extreme situations, alignment torque 
may also resTilt). (For a nonrotating star, the disk normal will be driven away from the stellar 
magnetic axis.) (§2.2). 

In general, we expect both types of torques to exist on the disk (§2.3 and §2.4). 

2.1. Diamagnetic Disk 

Here we consider the extreme situation where the disk is a perfect conductor and has no 
large-scale magnetic field of its own. The inner radius of the disk is located at r = r^, the 
magnetospheric radius. The magnetic field produced by the stellar dipole cannot penetrate the 
disk, and a diamagnetic surface current is induced. Aly (1980) has found the exact analytic 
solution to this model problem (see also Kundt Sz Robnik 1980; Riffert 1980). The magnetic field 
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at a point (r, (j),z = 0) (cylindrical coordinates) on the disk surface (z = 0, r > r^) is given by 

2/i 4/i 

Br = ^sinxcos((/)-0^ =F — tj^cosx, (2-1) 

^<t> = ^sinxsin((/)-(/>^), (2-2) 

= 0, (2-3) 



where x (/)^ are defined in Fig. 1 (they are both varying in time). In eq. ( |2-1| ), the upper 
(lower) sign applies to the upper (lower) disk surface, and the factor D is 

D = max (^^r^rl - 1, ^2H/r^ , (2-4) 

(where H <^ the half-thickness of the disk). The discontinuity in Bj- implies a surface current 

K = -^^cosx<A (2-5) 

{(j) is the unit vector along the (f) direction; similar notation will be used throughout the 
paper). Note that this surface current is induced to cancel the z-component of the stellar field, 
B^'^ = —iicosx/f^', the r, (/)-components of the stellar field [the first terms of eqs. ( |2-lj ) and ( |2-2| )] 
do not induce any net surface current. (In fact, they induce currents on the upper and lower 
surfaces of the disk, but these currents have opposite directions. Such currents can lead to the 
"squeezing" of the disk — changing the thickness of the disk, but not the "lifting".) The magnetic 
force per unit area on the disk results from the interaction between K and B^'^ (the first term in 
eq. 

^ = V¥7^^i^2xcos(</)-</.^)z. (2-6) 

The existence of this vertical magnetic force has already been noted by Aly (1980), and it is 
simply the difference in the magnetic pressure, B'^/{8tt), between the lower and upper surfaces. 
The magnetic torque per unit area is 

2u2 

N = rxF = J^sm2xcos{cj)-<P^)<p. (2-7) 

Clearly, averaging the force over the azimuthal angle, (F)<^ = (l/27r) Jq^ d(j)F, gives zero, but the 
azimuthally averaged torque is nonzero, and given by 

^2 2/^2 
^^)<^ = 2 ^ sm2x{sm (l)f,x -cos (l)f,y) = cos x (A x 0> (2-8) 

{fi and / are the unit vectors along the dipole moment and the disk angular momentum, 
respectively). For nonrotating star, this result implies that the magnetic torque tends to make the 
disk precess around the magnetic axis jl. For a rotating star — as long as the rotation period is 
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much shorter than the precession period, we need to average over the spin period. The fohowing 
identities wih be needed: 

cosx = cos /?cos — sin/?sin 6* sinwt, (2-9) 
sin X cos 0^ = sin 6* cos wt, (2-10) 
sin X sin 0^ = sin/3cos + cos /?sin sinwt, (2-11) 



where an appropriate phase for the stellar rotation has been adopted. We find that, after averaging 
over the spin period, the torque per unit 
Ps = 27r/a; is the spin period) is given by 



over the spin period, the torque per unit area ((N)) = ((N)<^)(^ = (l/Pg) Jq'" dt{N)^ (where 



2 

((N)) = ^^,^ cos/3f3cos^g- l) (lj X /). (2-12) 

The magnetic torque on a ring of radius r and width dr is simply dT = 27rr((N))dr. The angular 
momentum of the ring is (27rr$](ir)(r^O) (where S is the surface density of the disk, and 0, is the 
orbital angular frequency). Thus the effect of the magnetic torque is to make the ring precess 
around the spin axis cj at an angular frequency 

2 

fiprcc = ^2^7^^j^ COS (3 (3 COs2 9 - 1^ Cb. (2-13) 

Note that the sign of cos/3 (3cos^0 — 1) determines whether Clprec is along or opposite to it. 
This result has been obtained previously by Lipunov &; Shakura (1980) |^. 



2.2. Magnetically Threaded Disk 

We now consider the opposite limit in which the stellar magnetic field rapidly penetrates the 
disk (on a timescale shorter than the dynamical time of the disk) (see §2.3). Because of the shear 
between the disk and the plasma outside the disk, the threaded vertical field is winded to produce 
an azimuthal field which has different signs on the upper and lower surfaces of the disk. We thus 



^Lipunov & Shakura (1980) also argued that for cos^S < 1/3, the minimum of the interaction energy between the 
central dipole and the field generated by the disk current is achieved at f3 — 0° , while for cos^S > 1/3, the minimum 
energy corresponds to /3 = 90°. They therefore suggested that in the latter case (cos^^ > 1/3), the disk tends to 



evolve into the /3 = 90° state. However, if eq. (2-12) is the only torque present in the disk, it is not clear how /3 can 
change. Moreover, the "magnetic Bardeen-Petterson" effect always tends to align I and lu, independent of the sign 
of Qpiac (see §4.2). I became aware of the papers by Lipunov et al. (see also Lipunov, Semenov & Shakura 1981) in 
mid-February 1999, at which point this paper was nearly finished. For completeness and pedagogical reason, I have 
decided to keep this subsection (§2.1) in its original form. I thank Dr. Brad Hansen (CITA) for drawing my attention 
to the papers by Lipunov et al. 
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adopt the following ansatz for the magnetic field in the diskQ: 



Br = ^sinxcos((/> - 0^), (2-14) 
^<t> = ^sinxsin(0-</>^) ibC^cosx, (2-15) 
B, = -4cosx. (2-16) 



In ( 2-15| ), the second term represents the field produced by the twisting of Bz, and the upper 



(lower) sign corresponds to the value at the upper (lower) disk surfac^. The quantity C specifies 
the azimuthal pitch of the field line. In general we expect |C| ^ 1 (e.g., Sturrock & Barnes 1972; 
Lovelace et al. 1995 and references therein), but its actual value or form depends on details of the 
dissipative processes involved in the disk-magnetic field interactions. If the stellar magnetic field 
threads the disk in a closed configuration (e.g., Ghosh & Lamb 1979a,b; Wang 1987,1995), we 
expect C oc (O — w) so that C > ^or > a; and C < for 17 < w. But it has been argued that the 
differential shearing and the plasma flowing from the disk into the overlying magnetosphere will 
blow the field lines open and maintain them in a open configuration (e.g., Arons 1987; Newman 
et al. 1992; Lynden-Bell &: Boily 1994; Lovelace et al. 1995 and references therein), in which 
case we expect C, to be positive and of order unity. For our purpose in this paper, C, is simply a 
dimensionless number or function^. Also note that because of the screening currents in the disk 
(such as those discussed in §2.1) or in the magnetosphere boundary, the threaded field may be 
smaller than the vacuum field produced by the dipole. But since our result (see below) depends 
on we can easily absorb the screening effect into the definition of C- 



^ In principle, an additional radial field ABr could be generated (with opposite signs in the upper and lower 
disk surfaces) when the radial inflow of the disk drags the threaded B^. This discontinuous ABr will give rise to a 
precessional torque similar to that derived in §2.1. Since we expect ABr to be proportional to the radial velocity, it 



is neglected in eq. (2-14). 

^ There could be an additional contribution to B^ on the right-hand side of (2-15) due to the shearing of the 



stellar Br by the differential rotation. This contribution is not included for two reasons. First, the stellar Br may 
not be able to penetrate the thin disk; Second, even if a component of Br (e.g., the static component; see §2.3) 
penetrates the disk, the induced AB^ is expected to be smaller than that due to the threaded Bz. To see this, we 
write, schematically, 

^^^r(B.— +Br — ]-^, 
dt \ dz ^ dr ) Tdiss ' 

where Bz and Br are the threaded fields, Tdiss is the effective dissipation time. In steady state, we find 

3 

AB,^ — ^^(Bz — -Tdissf^Br, 

where ^ — {r / HB)Tdiss{^ — fj-") (for closed field configurations) or = {r / HB)Tdiss^ (for open field configurations) 
{Hb is the vertical scale in which Q varies). Note that the shearing of Br does not produce any surface current. For 
Hb ^ r, we can drop the term proportional Br in AB^, and obtain AB,^ — =f("Bz. 

^Even in models (e.g., Shu et al. 1994) where the disk is largely diamagnetic with no intrinsic Bz, the stellar field 
must penetrate the boundary layer. In this case, one would imagine that is nonzero only near the magnetosphere 
boundary. 
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The surface current on the disk corresponding to the field ( p-14D -( p-16 ) is 

Cc/i 



K 



cos X f. 



(2-17) 



The interaction of K with Bz gives rise to an azimuthal force, which tends to slow down (or speed 
up if c<j > and the field lines are in a closed configuration) the fluid motion, transferring angular 
momentum between the disk and the overlying magnetosphere or the star. This is the familiar 
magnetic breaking torque which has been included in previous studies of aligned magnetized disk 
(e.g., Ghosh & Lamb 1979a,b; Lovelace et al. 1995; Wang 1995; Lai 1998). Here we are interested 
in the vertical magnetic force (per unit area), resulting from the interaction between K and the 
azimuthal field produced by the stellar dipole (the first term in 2-15f| ): 



sin 2x sin( 



(2-18) 



This is simply the difference in B'^/ {8tt) below and above the disk. The magnetic torque associated 
with Fz is given by 



N 



sin 2x sin( 



Averaging over 



(N) 



we have 



sin 2x(cos (p^x + sin (l)^y) 



{i- fi) fi- {(L-i)i 



(2-19) 



(2-20) 



For a nonrotating star, this torque tends to pull the disk normal vector / toward being perpendicualr 
to the magnetic axis jl (assuming C > 0), thus making the disk plane parallel to the stellar field 
lines. For a rotating star, averaging over the spin period and using the indentities (|2^-( pll) , we 
find 



((N)) 



167rr^ 
87rr^ 



cos 



i2/3 (3 
/3(3 



cos 



1 



y 



cos 



{Co ■ 1)1 



(2-211 



In the absence of other forces, this magnetic torque will change the tilt angle /? of the ring (at 
radius r) according to 

Thus, depending on the sign of C cos/3(3 cos^ — 1), the angle /3 can increase or decrease: For 
C(3cos^0 — 1) > 0, the torque drives the disk toward the perpendicular configuration (/3 = 90°), 
while for ("(3cos^0 — 1) < 0, it drives the disk toward alignment (/? = 0°) or anti-alignment 
(/? = 180°). In particular, when C(3cos^0 — 1) > 0, the aligned configuration is unstable against 
the growth of disk tilt angle — This is the magnetic warping instability (see §2.3 and §2.4). 

It is instructive to understand the difference between the situation studied here and that of a 
diamagnetic disk considered in §2.1. In both cases, a vertical magnetic force is exerted on the disk. 
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Although eqs. ( p-6| ) and ( 2-18| ) appear similar, they have very different spin-averaged behavior. 
For a diamagnetic disk, we have 

2 

{F,)^ = sin 2(3 (3 cos2 6 - l) sin <p. (2-23) 

Thus the force on the y > side of the disk has a different sign from that on the y < side, giving 
rise to a torque which is along the x-axis (perpendicular to both I and cj). For a magnetically 
threaded disk, we have 

(F^)^ = ^^sin2/3(3cos2(9- 1) coscj). (2-24) 

The force has different signs for x > and x < 0, and the torque is along the y-axis (in the same 
plane as / and tD). 



2.3. A Hybrid "Realistic" Model 

The results of §2.1 and §2.2 represent two opposite, extreme situations, and are not likely to 
be realistic. The magnetic field configuration in a purely diamagnetic disk (see §2.1) is prone to 
Kelvin-Helmholtz instability and reconnection. It has been argued that the non-linear development 
of the these processes could lead to partial threading of the magnetic field through the disk (e.g., 
Ghosh & Lamb 1979a, b; Wang 1987, 1995), although the extent of the field-threading region is 
uncertain and the closed configuration advocated by Ghosh and Lamb is probably unrealistic 
(e.g., Arons 1987; Shu et al. 1994; Lovelace et al. 1995). In any event, near the magnetosphere 
boundary, the vertical stellar field must certainly thread the disk plasma. On the other hand, 
the magnetically thread disk model considered in §2.2 requires that the stellar field penetrates 
the disk almost instantaneously (compared to the spin period and the orbital period). This is 
unrealistic. The timescale for field penetration is uncertain, but it cannot be shorter than the 
dynamical time of the disk; it may be as long as the disk thermal time (E. T. Vishniac 1999, 
private communication; see Park & Vishniac 1996 and references therein). 

The vertical vacuum field produced by the stellar dipole on the disk can be written as the 
sum of a static component and a time-varying component: 

5^"^ = —-^ cos X = cos /3 cos 9 + sin/3sin0 sinwi. (2-25) 

While it is possible that the static field can penetrate the disk by allowing sufficient time for 
Kelvin-Helmholtz instability /reconnection to grow, it is almost certain that the variable field will 
be shielded out of the disk by the screening current. We therefore consider the following hybrid 
model, which we consider to be more realistic than the situations studied in §2.1 and §2.2: The 
static component of the vertical stellar field threads the disk, while the time-varying component 
is screened out by the disk. The winding of the threaded field will produce an azimuthal field 
and a radial surface current, while the variable vertical field will induce a shielding azimuthal 
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surface current and discontinuous radial field — this radial field can be obtained by appropriately 
modifying Aly's solution as given in §2. The magnetic field on the disk is then given by |^ 



Br = -J sinxcos((/> — 
= ■^sinxsin((/)- 
Bz = — cos /3cos 0. 



b^) ± (-^ cos P cos 



sin/?sin0sina;t, 



(2-26) 
(2-27) 
(2-28) 



In ( P-26|) , Bj. is the sum of the vacuum dipole field and the field produced by the azimuthal 
screening current (which is induced to cancel the variable part of Bi^^); in ( 2-27 ), B^ is the sum 
of the vacuum dipole field and the field induced by the winding of the constant B^ (see §2.2 for 
the property of C). The surface current on the disk is 



K 



2c/i 



sin P sin 9 sin ujt ( 



cos/? cos 6* f. 



The vertical magnetic force (per unit area) on the disk is then given by 



4/i2 



■ sin P sin 9 sin ut sin x cos( 



cos P cos 9 sin x sin( 



(2-29) 



(2-30) 



The corresponding torque (per unit area) acting on the disk, averaged over the ring, is given by 



2/i^ 



■ sin P sin 9 sin ut (fi x /) 



COS P COS 9 fi — {ji ■ 1)1 



(2-31) 



For nonrotating star, eq. ( 2-31| ) reduces to eq. (|2-20D (note that cosx = cos/? cos when lj = 0): 
The torque tends to make the disk plane align with the stellar field lines (i.e., / perpendicular to 
fi). Averaged over the stellar rotation, the torque can be written as 



((N)) = ((N))p,ec + ((N)Urp, 

where the precessional torque and the warping torque are given by 



((N))prec 
( (-^) ) warp 



cos/?sin^0 (w X /), 



.0^ 



sin 2P cos^9 y 



cos P cos^9 



UJ 



{Cj-l)l 



(2-32) 

(2-33) 
(2-34) 



Clearly, the diamagnetic feature of the disk induces a precessional torque (eq. | p-33| ), with the 
precession angular frequency given by 



prec 



cos P sin^0 UJ. 



(2-35) 



^In eq. (2-26) wc have neglected a possible component of Br generated by radial infall of the threaded Bz\ See 
footnote 2. Also, in (2-27) we have ignored a possible toroidal field generated by the shearing of the stellar Br; see 
footnote 3. 
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On the other hand, the disk-threading field gives rise to a torque (eq. [|2-34| ) which changes /? at a 
rate given by 

^ = sin 2/3 cos^^. (2-36) 

Note that the disk precession is along the direction (—cos /3a)), i.e., the precession is always in 
the opposite sense as the orbital motion of disk. Equation (|2-36| ) reveals the magnetic warping 
instability: the spin-orbit inclination angle /? increases when ^ cos /? > and decreases when 
Ccos/3 > 0. Thus, for C > (the most likely case; see §2.2), the warping torque always tends to 
drive the disk toward the configuration where I is perpendicular to cD. 



2.4. General Consideration 

The preceding subsections assume that the vertical field that threads the disk originates 
from the star. But it is important to note that this is not a requirement for the existence of the 
warping torque. Indeed, models of hydromagnetic driven outflows from disks are predicated on 
the existence of large-scale poloidal magnetic field that threads the disk — this field could have 
been advected inward by the accreting fiow or generated in the disk (see, e.g., Blandford 1989 for 
a review). 

Suppose there is vertical Bz (assumed to be time-independent) which threads the disk. The 
warping torque results from the interaction of the stellar i?^ and the radial surface current on the 
disk induced by the twisting of the threaded Bz- Thus in general, we can write the azimuthal field 
on the disk as 

^4: = ^ sin X sin((?i - 4>^J)T (Bz, (2-37) 

with ^ a positive dimensionless number of order or less than unity. The vertical force on the disk 
is given by 

Fz = ^smxsm{cP-cP^). (2-38) 
The averaged (warping) torque is then 

((N)) = ^^sin/3cos0y. (2-39) 

47]- 7-^ 

The model considered in §2.3 corresponds to Bz = — cos (3 cos 6 / r"^ . 

The precessional torque results from the interaction of the ring of diamagnetic screening 
current and the radial magnetic field produced by the star. Assume that the surface current has 
the form 

K = {K^i + K^2 sin tot) 0. (2-40) 

(In general, both the static current K^i and the time-dependent current K(f,2 smut are possible as 
the disk responds to the external field that tries to enters the disk.) The vertical magnetic force is 

Fz = - —^ {K^i + K^2smujt) sinxcos(<?!> - 0^), (2-41) 
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and the resulting averaged (precessional) torque is 

((N)) = (^ii:<^isin/?cos^ + ^cos/3sine^ x. (2-42) 

The model considered in §2.3 corresponds to Kipi = and K^2 = (2c/Li/7r^r^Z)) sin/3sin6', 
while in the model considered in §2.1, we have K^i = —{2c^/-jT^r^D) cos l3 cos 9 and 
K^2 = {2ctJ, / ir'^r^D) sin P sine. 

Finally, we note that in our calculations of the warping torques (§2.2-§2.4) we have assumed 
that the azimuthal pitch ( is stationary in time. An alternative scenario has been proposed (Aly & 
Kuijpers 1990; van Ballegooijen 1994), where the threaded field lines are being constantly wound 
up, with occasional field reconnection which releases the stored-up magnetic energy. Even for 
such a time-dependent magnetic field structure (apart form the time-dependence of the rotating 
dipole field), it is likely that magnetic torques (similar to those calculated in this paper) will exist, 
although this possible complication will be ignored in this paper. 



2.5. Disk Consisting of Diamagnetic Blobs 

The previous subsections treat the disk as a continuum. It has been suggested that under 
certain conditions, the accretion disk may be lumpy, consisting of diamagnetic blobs (Vietri & 
Stella 1998; see King 1993). Such inhomogeneous accretion flow may arise from the nonlinear 
development of various plasma instabilities (e.g., Kelvin-Helmholtz instability; Arons & Lea 1980) 
associated with the disk and the magnetosphere. Vietri & Stella (1998) studied the motion of 
diamagnetic blobs orbiting a central star with an inclined magnetic dipole and showed that 
magnetic drag force (Drell, Foley Sz Ruderman 1965) on the blob can induce vertical resonances 
near the corotation radius of the disk. 

We have calculated the magnetic torque on the accretion disk consisting of individual 
diamagnetic blobs. It can be shown that independent of the vertical resonances, there exists a 
magnetic torque which tends to induce tilt on the orbit of the blob. Since there are considerable 
uncertainties associated with such an extreme form of accretion disk, we relegate the calculation 
to Appendix C. 



3. MAGNETICALLY DRIVEN RESONANCES 

The discussion in §2 has neglected possible resonances in the interaction between the disk the 
rotating central dipole. Here we consider these resonances. 
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3.1. Vertical and Epicyclic Resonances 

The bending mode of the disk is characterized by a smah vertical displacement Z{r, 0, t), with 
the equation of motion 



where is the orbital angular frequency, the term —^^z represents vertical the gravitational 
restoring force (for a Keplerian disk, = and is the magnetic force (per unit area) as 
calculated in §2. Pressure, viscous force and self-gravity have been neglected. 

Consider first the magnetic force given by eq. ( |2-38 ) , and write it as 

^ sin6' COS sin(/) — (sin/?cos 6* -|- cos /?sin sinwt) cos (/> . (3-2) 



^ 27rr^ 



There exists the following possible vertical resonances: 

0^ = Q (for sin/?cos6' 7^ 0) (3-3) 

u:-VL = ±nz (for sin6'/0 and /? 7^ vr) (3-4) 

uj + Vt = ^z (for sin6'/0 and /? 7^ 0). (3-5) 

These resonances can be easily understood: (1) When sin /3 cos 7^ 0, there exists a static 
distribution (oc cos(/>) of vertical force field, and a fluid element feels the force oc cos fit; this 
explains the $7 = resonance. (2) When sm.6 7^ 0, there exists a time-dependent vertical force 
field (oc sinwt or coswt). A fluid element, traveling with angular velocity O, would experience 
vertical forces with frequency {lo it Vt) [for /? = 0, only (w — fi) is possible, while for (3 = 180°, only 
{(jj + VI) is possible]; this explains the a; it J7 = ztn^ resonances. 



Now consider the magnetic force given by eq. ( 2-41j ). The force associated with K^i has the 



form sm.(j) and cos((/) it wt), which would give rise to the Vl^ = Vt resonance and the a; ± = itri^ 
resonances as in eqs. ( p-3| )-( |3^ ). The force associated with K^2 bas a time-dependence of the form 



sinwt sinxcos((/> — (/)^) = — cos /3 sin sin (/) -|- sin /? cos sin sin < 



2 

This gives rise to the following resonances: 



H — sin 9 {sm. 2u)t cos (j) — cos (3 cos 2u!t sin (f)^ . (3-6) 



nz = ^ (for cos /3 sine 7^ 0) (3-7) 

uj±n = ±nz (for sin /3 cos 6* 7^ 0) (3-8) 

2uj-n = ±nz (for sin 6* 7^0 and /? 7^ vr) (3-9) 

2uj + n = nz (for sin6l7^0 and /? 7^ 0). (3-10) 



The new resonances, 2lo ziz Q = ±0^, come about because the magnetic field varies as costot or 
sinwt, and the screening current also varies as cos ujt or sinwt. 
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In addition to the vertical resonances discussed above, epicyclic resonance can also arise 
from the magnetic field - disk interaction. Consider the field given by eqs. ( |2-26| )-( |2-28| ). The 
interaction between with gives rise to a radial force 

Fr = - 2^2^6^ sin 2/3 sin 26* sin a;t. (3-11) 

Clearly, epicyclic resonance occurs when 

K = uj (for sin 2/3 sin 26* / 0), (3-12) 

where k is the epicyclic frequency. For a Keplerian disk, k = Q, ( |3-12D is a corotation resonance. 

We note that this paper deals only with dipole field from the star. When higher-order 
multipole fields are considered, it is conceivable that additional resonances can arise. 



3.2. Magnetic Torques at Resonances 



What are the consequences of the magnetically driven resonances? We have not investigated 
this issue in detail. The resonances may act as an extra source (in addition to the non-resonant 
warping torques discussed in §2 and §4) for generating bending waves and spiral waves in the disk. 
Near the resonances, fluid elements undergo large out-of-plane and radial excursions, which may 
lead to thickening of the disk^. 

We can get some insight into the resonances by calculating the magnetic torques at the 
resonant radii. Consider first the torque associated with the force in eq. (3-2). For a Keplerian 
disk, the condition = is always satisfied. When sin /3 cos ^ 0, there always exists a net 
(averaged) torque acting on the fluid element, pulling its orbital angular momentum axis toward 
(or away from) the spin axis (eq. ||2-39(| ). The condition + fi^ = w (or 20 = u for a Keplerian 
disk) specifies a special radius in the disk. To calculate the net torque on the fluid element at this 
resonant radius, we cannot average (p and the spin period independently (as done in §2). Instead 
we write the torque N = r x F as: 



N 



47rr^ 



sin 9 cos ujt{l — cos 2<j)) — (sin (3 cos 9 + cos /? sin 9 sin cot) sin 20 



+ 



sin 9 cos Lot sin 2(j) + (sin /3 cos 9 + cos /? sin 9 sin ut){l + cos 2</)) y 



(3-13) 



Following a fluid element we have 
resonance, the torque is given by 



((N)) 



sin (3 cos 9 y 



^It + (pQ. Averaging over time, we find that, at the 20 = uj 



(1 + cos /?) sin 9 {x cos IcpQ + y sin 2(/)o 



(3-14) 



®This may be analogous to the Lorentz resonances (which occur when charge particles move around a rotating 
magnetic field) in the jovian ring (Burns et al. 1985; Schaffer & Burns 1992). However, because of the fluid nature 
of the disk, the resonances may not lead to sharp edges in the disk. 
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Thus for each fluid element on the resonant radius, the averaged torque is modified from the 
"nonresonant" value (the first term in eq. p-14| ). Since different fluid elements in the same ring 
have different values of (pQ, they will experience different torques; if they are allowed to move 
independent of each other, the ring will eventually disperse. However, if there is strong coupling 
between the different elements on the ring so that the ring evolve dynamically as an identity, 
then one should average over </)o, and eq. ( 3-14| ) reduces to (|2-39|) . We can similarly consider the 
torque associated with the magnetic force given in eq. ( 2-4l| ). At both the = uj and 217 = uj 
resonances, in addition to eq. ( ^-42 ), there are "resonant" torques which depend on (pQ (analogous 
to eq. [FTP - 

We note that the "nonresonant" torque is nonzero only when sin/? ^ (misaligned spin-orbit), 
while at the 20 = u; or = a; resonances, the torque can be nonzero even when /3 = 0. Also note 
that in general, the "resonant" torque is of the same order of magnitude as the "nonresonant" 
torque. We shall not consider these resonances in the rest of the paper. 



4. DYNAMICS OF WARPING AND PRECESSING DISK 

We now study the dynamics of the disk under the influence of the magnetic torques calculated 
in §2. For concreteness, we shall use the torque expressions of §2.3; using other expressions of 
§2 would give similar results (although the dependence on angles would be different). We are 
particularly interested in whether the warping instability can operate in the presence of disk 
viscosity. 



4.1. Criterion for the Warping Instability 



Our starting point is the evolution equation for disk tilt l{r,t) (the unit vector perpendicular 
to the disk annulus at radius r) given by Pringle (1992) (see also Papaloizou & Pringle 1983) [|: 



dl 



Vr 



d (i 



dl 



1 



m 

dr 



m 

dr 

2 



+ 



N 



(4-1) 



Here Vr is the radial velocity of the flow, 17' = dQ/dr, ui is the usual disk viscosity (measuring the 
r-(j) stress), and V2 is the viscosity (measuring the r-z stress) associated with reducing disk tilt. 
We assume that the timescale for I to change is much longer than the spin period, so that N is the 



'^Recent study (Ogilvie 1999) indicates that in the nonhnear regime, the warping equation needs to be modified. 
In effect, vi and V2 are not the usual vertically averaged viscosities and may depend on the amplitude of the warp. 
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averaged torque (per unit area) as calculated in §2.3 (The notation ((• • •)) has been suppressed). 
Using the following relations for a Keplerian disk: 



Vr 



' 2r ^ ' 



M 



J, 



(4-2) 



where ^ is a function of r which approaches unity in the region far from the inner edge of the disk 
(see Appendix A), and assuming that 1^2/^1 is independent of r, we reduce eq. ( [4-l| ) to 



m 

dt 



2rJ'\ 3i^i 



1 d^l 1 



+ 



2r 



J 



-1 



1 



m 

dr 



dl 



dr 



+ 



N 



(4-3) 



We rewrite the magnetic torques in eqs. ( 2-33| ), ( p-34 ) as 
N 



" prcc 
Nwarp 



-Qp cos P Lj X I, i7 
cos (3 sin (3 y, T 



= ^ 



sm 



cos 



(4-4) 
(4-5) 



At this stage / is no longer useful as a coordinate axis, so we consider a different cartesian 
coordinate system where to is the z'-axis (see Fig. 1). In this coordinate system, we write 

/ = (sin /3 cos 7, sin fj sin 7, cos /?) , (4-6) 

which defines the local twist angle 7. For small tilt angle (/3 ^ 1), eq. ( |4-3| ) simplifies to 

dW 3u2dW 1 d'^W 



dt 



4r c?r 2 ^ Sr^ 



(4-7) 



where W = /?e*'^, and we have set = 1 (far from the inner edge of the disk) - using a more 
rigorous J would only affect the numerical coefficient in front of the dW/dr term. 

To derive the instability criterion, we consider the WKB solution of the form (valid for 
kr > 1) 

W (X exp{iat + ikr) . (4-8) 



Substituting into 



we obtain the dispersion relation 



a 



3l^2 

4r 



k -Qp \ + i[ \i'2k^ - ^ 



For the instability to grow, we require Im((T) < 0, or 

1 



Fu, > 2^2^^ *^=^ Instability. 



(4-9) 



(4-10) 
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Since the radial wavelength is restricted to A < r, or A; > 27: /r, the instability criterion becomes 



r^„ > 27r2^ Instability. (4-11) 



This criterion has a simple physical interpretation: The magnetic torque drives the growth of disk 
tilt on a timescale T"^, while viscosity tries to reduce the tilt on a timescale jv^ (which is of the 
same order as the radial drift time of the flow, r/|I4-|) if is of the same order as v\). Instability 
requires T"-"^ <, r"^ /v2- 



Using (4-2) and ( [4-5| ) , and assuming that 1^2/1^1 is independent of r, we can further reduce 

(Eni) to 

This indicates that inside the critical warping radius r^, the magnetic torque can overcome the 
viscous force and make the disk tilt grow. Moreover, with the expected^ C ~ 1; v\l^2 ~ 1 and 
i7 < 1, the warping radius is of the same order of magnitude as the magnetosphere radius (see 



eq. [1-1] and Appendix B for a discussion of the magnetosphere radius for arbitrary geometry). 
Thus, as the disk approaches the magnetosphere, it will tend to be tilted with respect to the 
stellar spin even if at large radii the disk normal is aligned with the spin axis. 

Similar analysis can be applied to the case where the disk is nearly antiparallel to the spin 
(vr — /? <C 1). If we define W = {tt — P)e^^ , then eq. ( [4-71 ) remains valid except the sign in front of 
ii}pW changes. We obtain the same instability criterion as above. Thus, an antiparallel disk tends 
to be driven toward a perpendicular configuration near the magnetosphere radius. 



4.2. Magnetic Bardeen-Petterson Effect 

A tilted disk will be driven into precession by the torque Nprec with fiprec = —^pcos (3 lj (see 
eq. [}4-4|| ). What is the effect of this precession on the disk tilt? An analogy can be made with 
the behavior of a disk undergoing Lense-Thirring precession around a rotating black hole (or any 
compact object) (Bardeen & Petterson 1975; Kumar & Pringle 1985; Pringle 1992; Scheuer & 
Feiler 1996). 

The gravitomagnetic force from the rotating object (with spin angular momentum J = Jui) 
drives the precession of the misaligned disk at angular frequency f^LT = 2GJ/ (c^r^). Bardeen 
& Petterson (1975) pointed out that the action of viscosity on the differentially precessing 
disk tends to align the rotation of the inner disk with the spin axis of the central object. The 



*For Keplerian and inviscid (or very nearly so) disks, resonance between the epicyclic frequency and orbital 
frequency leads to 1/2/ui — l/(2o?) ^ 1. However, the resonance is delicate and might be destroyed by the effects of 
general relativity, self-gravity, magnetic fields and turbulence. See Ogilvie (1999) for a discussion. 



- 18 - 



Bardeen-Petterson radius tbp, that is, loosely speaking, the radius inside which the disk aligns 
with the spin, is obtain by equating the precessional time J7£"rJ and the viscous time r'^ /v2 (Scheuer 
& Feiler 1996), i.e., 

TBP = (4-13) 

However, the transition from the warped outer disk to the aligned inner disk is rather broad, and 
for this reason the timescale to achieve the Bardeen-Petterson alignment is much larger than the 
precession timescale at rep (Kumar & Pringle 1985; Pringle 1992). Numerical simulations (Pringle 
1992) indicate that a steady state is achieved on a time scale of order (10 — 100) f^LT (evaluated 
at TBp). 

Now consider the effect of magnetically driven precession^. Setting 17^ equal to 1^2 /r'^-, we 
obtain the magnetic Bardeen-Petterson radius: 

Inside tmbP; the combined effect of viscosity and precession tends to align the disk normal with 
the spin axis. We see that typically tmbp is of the same order as (the warping radius) and 
(the magnetosphere radius). Thus the precessional torque has an opposite effect on the disk tilt 
as the warping torque discussed in §5.1. However, because of the broad warp-alignment transition 
expected for the magnetic Bardeen-Petterson effect and the long timescale involved, we expect 
that the precession-induced alignment will be overwhelmed by the warping instability. 



5. EFFECT ON THE SPIN EVOLUTION 

How does the magnetically warped and precessing disk affect the spin evolution of the 
central star? The angular momentum of the accreting gas is deposited at the magnetospheric 
boundary and transferred to the the central object. In addition, there are back-reactions on the 
star associated with the magnetic torques on the disk. Thus the spin angular velocity u evolves 
according to 

^ {lu) = M{GMr^)^'n - A/'prec - A/'warp, (5-1) 

where / is the moment of inertia, A/'prec and A/'warp are the total warping and precessional torques 
acting on the disk |^. To calculate A/'prec and A/'warp i we need to know the tilt angle /3 as a 



'in some systems (such as accreting neutron stars with weak magnetic fields), the Lense-Thirring precession 
dominates (See §6). But here we neglect SIlt- 



^ In (5-1 ) we have neglected possible magnetic torque [as in the Ghosh & Lamb (1979a, b) picture] associated closed 
magnetic field lines that thread the disk and connect the star. Such magnetic torque can be included by modifying 
the first term on the right-hand-side of (5-1) to M^JGMrm f I, where / is a dimensionless function which depends 
on the ratio uj/Q{rm) [See, e.g., eq. (35) of Lai (1998) and references therein]. 
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function of r. In principle this can be obtained by solving the tilt equation (4-1) or ( |4-3| ) . But note 
that since the torques per unit area (eqs. [ p-33 |-| 2^^3^ ) are steep functions of r, we can assume, as 
a first approximation, that [3 is independent of r near r^- We then find 



and 



A/", 



warp 



prec 



27rrNwarp dr 



27rrNprec dr 



sm 



'cos/?(a) X /), 



■ COS a sm 2p y 



I2rl 



■ —r- COS (J COS p 



(cj • 1)1 



(5-2) 



(5-3) 



Note that the ratio of the characteristic magnetic torque, A^mag = /^^/'^m; ^-^d the characteristic 
accretion torque, A^acc = M^/GMr^, is 



mag 



7/2 



V 



'7/2 



[see eq. (1^) and Appendix B]. The spin evolution equation can then be written as 

, CCOS^0 . 2r, n - 

•I ^ acc 

rccos^6' 2^ J . ^ . 

+A/acc g^yy^ cos - 1 sm/? u;_L 



(5-4) 



cos [3 - 

2^cos^ ) / X (2), 



sm 



(5-5) 



V37r?77/2 

where u)^ is a unit vector perpendicular to a) and lies in the plane spaned by / and a) (see Fig. 1 



The physical effects of the three terms on the right-hand-side of (5-5) are evident: The first term 
is responsible to the spin-up or spin-down (depending on the sign of cos /3) of the star; the second 
term tends to align or misalign the spin axis with the disk axig|^, and the third term induces 
precession of the star's spin axis uj around /. The timescales associated with these changes of u 
(spin-up/spin-down, alignment/misalignment and precession) are all of order 



' spin 



A/" ICC 



luj 



M./GMF„ 



(5-6) 



This spin-changing timescale is typically much longer than the timescale associated with changing 
the disk orientation (see §4 and §6). 

Note that eq. (|5-6| ) represents the instantaneous spin- up/spin-down timescale (i.e., at a given 
f3). If the disk inclination /? wanders around 90°, we expect the secular spin-up of the star to be 
slower (see §6). 



Since the timescale to change ili is much longer than the timescale to change I, the change of the tilt angle /3 is 



determined by the dynamics of / rather than lo (see eq. [2-36| and §4) 
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6. ASTROPHYSICAL APPLICATIONS 

In this section we discuss/speculate several possible applications of our theory. While 
preliminary, they demonstrate the potential importance of the physical effects uncovered in 
previous sections. 

6.1. Spin Evolution of Disk- Fed X-Ray Pulsars 

Recent long-term, continuous monitoring of accreting X-ray pulsars with the BATSE 
instrument on the Compton Gamma Ray Observatory has revealed a number of puzzling behaviors 
of the spins of these accreting, magnetized neutron stars (see Bildsten et al. 1997 and references 
therein). Several well-measured disk- fed systems (e.g., Cen X-3, GX 1+4: and 4U 1626-67) display 
sudden transitions between episodes of steady spin-up and spin-down, with the absolute values of 
spin torques approximately equal (to within a factor of two). Of special interest is the observed 
anticorrelation between the torque and X-ray luminosity during the spin-down phase of GX 1-1-4 
(i.e., the torque becomes more negative as the luminosity increases; Chakrabarty et al. 1997). 
These features are at odds with previous theoretical models, according to which the neutron star 
must be near spin-equilibrium in order to experience both spin-up and spin-down (e.g., Ghosh &: 
Lamb 1979a,b; see also Yi, Wheeler &; Vishniac 1997; Torkelsson 1998; Lovelace et al. 1998). It 
has been noted that the observational data can be nicely explained if the disk can somehow change 
its sense of rotation (Nelson et al. 1997) — This poses a significant theoretical problem, since the 
disk formed in a Roche lobe overflow has a well-defined direction of rotation. It has been suggested 
(Van Kerkwijk et al. 1998) that the disk reversal may be caused by the radiation-driven warping 
instability (Pringle 1996), although the extent and timescale of the reversal remain unclear. 
(Recall that one expects the radiation-driven instability to operate only at large radii; Pringle 
1996.) Also, since the radiation driven warping does not directly depends on the spin orientation, 
it is not clear why the disk prefers to wander around the perpendicular state (/3 = 90°). 

Here we suggest that the magnetically driven warping instability uncovered in this paper 
plays an important role in the determining the spin behaviors of accreting X-ray pulsars. The 
magnetosphere is located at 

/ 4 \ V7 

and we shall consider 77 to be a constant of order unity (see Appendix B). Here /lisq, Mi 4 and M17 
are the neutron star's magnetic moment, mass and accretion rate in units of 10^° Gcm^, 1.4M0 
and 10^^ g s^^, respectively. Near the magnetosphere, the disk lies in the so-called "middle" 
(gas pressure and scattering dominated) region of the a-disk solution (Shakura &; Sunyaev 1973; 
Novikov & Thorne 1973). The surface density of the disk is 

S = (1.7 X 103 g cm-2) «-4/5M;//Mff r8-'/V'/^ (6-2) 
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where rg = r/(10^ cm), and the possible effect of a threaded magnetic field can be included in the 
definition of J (see Appendix A). The growth rate of the magnetic warping instability (eq. [ [4-5| ]) 
is given by 

= (0.035 s-i) (C cos2 0) aVs^l^MrrVj-'/Vs^'/'V-^/^ 

As discussed in §2 and §4, the magnetic torque tends to drive the inner region (near the 
magnetosphere) of the disk toward a perpendicular configuration (/? = 90°) on a timescale of order 
r~^. Indeed, the (3 = tt/2 state represents an "attractor". In the idealized situation (each disk 
ring evolves independent of each other), the tilt of a nearly orthogonal disk evolves according to 

^ = sin /? cos /3 ~ -r^ A/3, (6-4) 

where A/3 = /? — t: /2 (see eq. |]2-36| ). The star can then spin-up or spin-down, depending on the 
whether /3 < 90° or /3 > 90° [see eq. (|5^)]: 



1 H TTj:— sm (3 



cos/3. (6-5) 



6^7/2 

The characteristic spin- up/spin-down time scale (eq. [ p-6( |) is 

r,pi, = (7.9 X 103 yrs) r'''M-fh, fi^.'^' M-'/' (i^) , (6-6) 

where Ps is the spin period, and = //(lO^^ g cm^). 

In our picture, the observed sign switching of lo in several X-ray pulsars is associated with 
the "wandering" of (3 around the "preferred" value (/3 = 90°). Such "wandering" needs not be 
periodic: Consider a disk initially at /3 = 90°. Imagine that a perturbation in the accretion 
induces a negative (not necessarily small) A/3 = (3 — 90° — This is achieved on a viscous timescale 
(which is of the same order of magnitude as F"^ at the inner disk edge; see §4.1). The star 
spins up. The magnetic torque then drives (3 toward 90° on the timescale of T"^, at which point 
another perturbation (which cannot be faster than T"^) can induce another A/3, which can be 
either negative or positive, and the star will then continue to spin- up or switch to spin-down. 
We expect that the timescale of the switching between spin-up and spin-down is of order a few 
times (evaluated at the inner disk boundary). For a ~ 0.1 — 0.01 and Mn ~ 0.1, /iso ~ 1 
(typical of X-ray pulsars), eq. ( |6-3| ) gives F~^ ~ 5 — 30 days (for Ccos^O = 1), comparable to what 
is observed in Cen X-3 {Ps = 4.8 s) (Bildsten et al. 1997). For 4U 1626-67 {Ps = 7.6) and GX 
l-|-4 (Ps = 120 s), the sign of uj switchs once in 10 — 20 yr. Since equation ( |6-3[ ) is uncertain and 
depends on many parameters, it is conceivable that such long switching time can be accommodated 
(e.g., with Qcos^O = 0.1, a = 0.01, M17 = 0.1 and /U30 = 10, we find F~^ ~ 6 years). 
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Note that eq. ( |6-6| ) represents the "instantaneous" spin-up/spin-down time of the neutron 
star. Since the outer disk has a well-defined direction (presumably with f3 < 90°), the inner disk 
will unlikely spend equal amount of time in the prograde (/? < 90°) phase and the retrograde 
(/3 > 90°) phase. Thus we expect that on a longer time scale (longer than the disk reversal 
time r~^), the neutron star will experience secular spin-up. This expectation is borne out by 
observations, although it is difficult to predict the long-term spin-up rate. 

A full study or simulation of the nonlinear behavior of the inner disk tilt would be desirable 
to make more meaningful comparison with observations. But our discussion and estimate given 
above indicates that magnetically driven warping may be a crucial ingredient in explaining the 
spin behaviors of disk-fed X-ray pulsars. Other physical effects (such as radiation driven warping 
and propeller effect; Van Kerkwijk et al. 1998, Lovelace et al. 1998) may also play a role. 

6.2. Weakly Magnetized Neutron Stars: Quasi-Periodic Oscillations 

Rapid variability in low-mass X-ray binaries (LMXBs), containing weakly magnetized 
(B ~ 10^ — 10^ G) neutron stars, has been studied since the discovery of the so-called horizontal- 
branch oscillations (HBOs) in a subclass of LMXBs called Z sources (van der Klis et al. 1985; 
Hasinger & van der Klis 1989). The HBOs are quasi-periodic oscillations (QPOs) (with the 
Q- value v j of order a few, and rms amplitude ^ 10%) which manifest as broad Lorentzian 
peaks in the X-ray power spectra with centroid frequencies in the range of 15 — 60 Hz which are 
positively correlated with the inferred mass accretion rate (see van der Klis 1995 for a review). 
For many years, the standard interpretation for the HBOs has been based on the magnetosphere 
beat-frequency model, first advocated by Alpar & Shaham (1985) and Lamb et al. (1985), in 
which the HBO is identified with the difference frequency between the Keplerian frequency at the 
magnetospheric boundary and the spin frequency of the neutron star (see Ghosh & Lamb 1992 
for a review). However, recent observations of kHz QPOs (500 — 1200 Hz) in at least eighteen 
LMXBs by the Rossi X-ray Timing Explorer (RXTE) have called into question this interpretation 
of HBOs (see van der Klis 1998a, b for review). The kHz QPOs (with Q up to 100 and rms 
amplitude up to 20%) often come in pairs, and both frequencies move up and down as a function 
of photon count rate, with the separation frequency roughly constant (The clear exceptions are 
Sco X-1 and 4U 1608-52; van der Klis et al. 1997, Mendez et al. 1998a). In most Z sources, the 
15-60 Hz HBOs appear simultaneously with the kHz QPOs, while in several atoll sources (which 
are thought to have weaker magnetic fields and smaller accretion rates than the Z sources), broad 
peaks at 10 — 50 Hz in the power spectra (similar to the HBOs in the Z sources) have also been 
detected at the same time when the kHz QPOs appear. While the origin of the kHz QPOs is 
uncertain, it is natural to associate the higher frequency QPO with the orbital motion at the 
inner edge (perhaps the magnetosphere boundary) of the accretion disk, and the lower-frequency 
QPO may result from the (perhaps imperfect) beat between the Kepler frequency and the neutron 
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star spin frequency]^ — This beat frequency interpretation is supported by the observations of a 
third, nearly coherent QPO which occurs during X-ray bursts, at a frequency approximately equal 
to the frequency difference between the twin kHz peaks or twice that value. (An exception is 4U 
1636-53, Mendez et al. 1998b.) Clearly, if this generic identification of kHz QPOs is correct, the 
beat between the ~ 1 KHz Keplerian frequency at the disk inner edge and the ~ 300 Hz spin 
frequency cannot produce the 10 — 60 Hz low-frequency QPOs (LFQPOs: HBOs in the Z sources 
and similar features in the atoll sources), unless one postulates (Miller, Lamb & Psaltis 1998) that 
the inner disk edge lies inside the magnetosphere and is unaffected by the magnetic field. 

Stella and Vietri (1998) suggested that the 10 — 60 Hz LFQPOs are associated with 
Lense-Thirring precession of the inner accretion disk around the rotating neutron star. Assuming 
that the low-frequency QPO and the kHz QPOs are generated at the same radius in the disk, 
one obtains, to leading order, Qlt = (2/a;/Mc^)i7^ (where Q,, Qlt and to are the orbital, 
Lense-Thirring, and spin angular frequencies, respectively). Thus the LFQPO frequency depends 
quadratically on the kHz QPO frequency. (The classical precession associated with spin-induced 
oblateness of the star, as well as relativistic effect beyond the Lense-Thirring formula, both can 
introduce small correction to this correlation; see Stella &: Vietri 1998, Morsink & Stella 1999.) 
This is consistent with observations of several sources (e.g., the Z sources GX17-I-2, GX5-1, 
Sco X-1, and the atoll source 4U1728-34; see Ford & van der Khs 1998, Psaltis et al. 1999 and 
references therein) P^ . 

For the Lense-Thirring interpretation of HBOs to be viable, the inner disk must be tilted 
with respect to the stellar spin axis. The Bardeen-Petterson effect tends to keep the inner region 
of the disk [inside tbp, typically at (100 — 1000)GM/c^] co-planar with the star (Bardeen & 
Petterson 1975). Radiation driven warping (Pringle 1996) is only effective at large disk radii. 
While global disk warping modes may exist with nonzero tilt near the inner disk boundary (Ipser 
1996; Markovic & Lamb 1998), an external driving force is needed to excite them. 

Here we suggest that the magnetic warping torque provides a natural driver for the disk tilt 
near the inner accretion disk |^. For typical parameters of LMXBs, the magnetosphere is located 



^^See Stella & Vietri (1999) for an alternative interpretation of the lower QPO peak which does not involve beating. 

^''Observations indicate that the ratio I/M required to fit the expected Olt — ^ relation is a factor of 2 — 4 larger 
than allowed by neutron star equation of state. This situation can be improved if the observed LFQPO frequency is 
the second harmonic of the fundamental precession frequency (See, e.g., Stella & Vietri 1998, Morsink & Stella 1999, 
Psaltis et al. 1999). In some Z sources, one requires that the observed HBO frequency is four times the fundamental 
precession frequency in order to produce reasonable I/M. Alternatively, the spin frequency is twice of what is inferred 
from the difference between the twin kHz QPOs — This would make the beat frequency interpretation of the lower 
kHz peak invalid. It would be interesting to search for "sub-harmonic" feature of the HBOs in the power spectra (see 
Ford & van der Klis 1998 for possible evidence of such a sub-harmonic feature in 4U 1728-34). 

^■^Vietri and Stella (1998) suggested that if the accretion disk is inhomogeneous, diamagnetic blobs can be lifted 
above the equatorial plane through resonant interaction with the star's magnetic field near the the corotation radius. 
See Appendix C. 
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at0 

Tm = (18 km) r, M^l" M^,^/' (6-7) 

where /U26 = /^/(lO^^ G cm^). Assuming that near the magnetosphere the disk is described by the 
"inner region" (radiation and scattering dominated) solution of Q-disk (Shakura & Sunyaev 1973; 
Novikov & Thorne 1973), we have for the disk surface density and half-thickness 

S = (l05 g cm-2) a-iMi7]/^Mf7V^/V"\ (6-8) 
H = (l.lkm)Mi7J, (6-9) 

where rg = r/(10^ cm). The growth rate of the magnetic warping instability (eq. ]4-5[ ) is then 

= (555s-^) (Ccos^e) a/x^6^i7rg~^:7 

1.0 s-i) (Ccos^e) M,.,^,-,iMl, (^)'^- ^^-^^^ 



The magnetic effect also results in a precessional torque on the disk. From eq. ( 2-35| ) and using 



8[)-( |6-9| ), we find that the precession frequency associated with this magnetic torque is 
i.p,ec = - (0.21 Hz) (sin^e cos /?) M1.4 n^iufj ' D'^J 



(0.60Hz) (sin^^cos/?) (^^] mI'J mI^'^^^ (^)'^'^'' ^^-11 



where in the second equality we have used D = {2H /vmY^'^ [see eq. ( ^-4] )]. This should be 
compared with the Lense-Thirring precession frequency 



i/LT = (44.4Hz) /45r, 



-3 



300 Hz 



= (8.1 Hz) Ml'2 /45 ,~,r''Mlp (^) (^) , (6-12) 
where is the spin frequency. Note that for [3 < 90° the magnetically driven precession is opposite 



to the stellar spin axis (thus the negative sign in eq. [ 5-ll| ), while Lense-Thirring precession is in 
the same direction as the spin. For a given source, the magnitude of |fprec| is typically smaller 
than z/LT, but it becomes increasingly important with increasing M (since z^lt oc M^/"^ while 
i^prec oc M^'^/^). This may explain the observed flattening of the correlation between the LFQPO 
frequency and kHz QPO frequency as the latter increases |^. 
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Note that since Vm is so close to the inner-most stable orbit, general relativistic effect tends move the inner disk 



edge to a radius larger than what is given in (6-7); see Lai (1998). Here we shall neglect such complication. 

^^Since the classical precession rate due to the oblateness of the star is negative (for cos /3 > 0) and scales as r~''l'^ 
(while i^LT oc r~^) (Stella & Vietri 1998, Morsink & Stella 1999), it may also explain this observed flattening, provided 
that the spin frequency is much higher than inferred from the kHz QPOs and the burst QPOs. However, because 
of the similar power-law indices (r~'' ^"^ vs. r~^), it is difficult to explain why the !/lfqpo oc z'kHzQPO scaling breaks 
down only at very high kHz QPO frequency (Psaltis et al. 1999). 
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We note that although our quahtative conclusion that magnetic effect can induce tilt of the 
inner accretion disk and therefore facilitate its precession is robust, the analytical expressions 
given in this subsection only serve as an order-of-magnitude estimate which indicates the potential 
importance of the magnetically driven precession. More detailed calculations (including global 
mode analysis) are needed to make more meaningful comparison with observational data. In 
addition to the intrinsic uncertainties associated with the a-disk model and the magnetosphere 
boundary layer, there are two complications that may prove important for LMXBs: (i) General 
relativistic effect can modify the inner radius of the disk so that it is not just determined by 
magnetic-plasma stress balance (Lai 1998); (ii) The magnetic field is not expected to be dipolar. 
This comes about either because of the multipole fields from the central star, or, even if the 
intrinsic stellar field is dipolar, the (partially diamagnetic) disk can enhance the field in the 
boundary layer (Lai, Lovelace & Wasserman 1999). 



6.3. Super-Orbital Periods in X-ray Binaries 

Here and in §6.4 we speculate upon two additional applications of our theory. 

Many X-ray binaries are known to exhibit long-term cycles in their X-ray or optical 
luminosities (see Priedhorsky &; Holt 1987 and references therein). Of particular interest is the 
well-established "third" period (longer than the orbital period) in Her X-1 (35 d) (Tananbaum et 
al. 1972), LMC X-4 (30.5 d) (Lang et al. 1981) and SS433 (164 d) (Margon 1984). It is generally 
thought that these super-orbital periods result from the precession of a tilted accretion disk. For 
Her X-1, Katz (1973) proposed that the precession was forced by the torque from the companion 
star, but left unexplained the origin of the disk's misalignment with respect to the orbital plane. 
Recently it has been suggested that the radiation-driven warping instability (Pringle 1996) is 
responsible for producing warped, precessing disks in many X-ray binaries which exhibit long-term 
cycles (Maloney &: Begelman 1997; Wijers & Pringle 1998). 

It is likely that magnetic torque plays a role in driving the warping of the inner disk. The 
observed systematic variation of the X-ray pulse profile of Her X-1 requires the inner edge of the 
disk to be significantly warped (Sheffer et al. 1992; Deeter et al. 1998 and references therein). In 
addition, the magnetically driven precession frequency is 

where is the magnetosphere radius [eq. (|6-1| )], and we have adopted the "middle region" 
solution of the a-disk (near the inner edge of the disk, D ~ 0.2). This is comparable to the 
observed super-orbital periods. Moreover, the magnetically driven precession is retrograde with 
respect to the direction of rotation of the disk, in agreement with observations. Of course, the 
precession rate is a function of r and /?, so a modal analysis is needed to determine the global 
precession period. 
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6.4. T Tauri Stars 

It is well established that classical T Tauri stars (CTTS) have circumstellar disks; Evidence 
also exists for magnetospheric accretion induced by the stellar magnetic field (e.g., Hartmann 
1998). For typical parameters (M ~ 0.5 M©, R ~ 2Rq, M = - W'^ Mq yr~\ and surface 
field ^ 1 kG), the magnetosphere is located at a few stellar radii. Using n = B^^R^, we find 
that the growth time for the warping instability at the inner region of the disk is given by [see 
eq. dU)] 

r- = (6.5day.)(— ) (^) (— ) (— ) (^__,) (ccos^ «) . (6-14) 



The precession period of the tilted disk (see eq. | 2-35|] ) is of the same order of magnitude as F"^. 
The surface density S is unknown, but reasonable estimates give Tj 1 — 100 g cm^^ (Hartmann 
1998). Therefore F~^ ranges from days to years, much shorter than the lifetime of T Tauris stars 
(~ 10^ years). 

It has been observed that the photometric periods (5 — 10 days) of CTTS vary (by as much as 
30%) on a timescale of weeks (Bouvier et al. 1995). The origin of this variability is unknown. If we 
interprate the photometric period as the orbital period at the magnetosphere boundary (Bouvier 
et al. 1995), then we may understand the period variation in the context of warped, precessing 
disks: As the inner disk warps and precesses, it experiences different stellar magnetic field, and 
thus the inner disk radius varies. 

It is also of interest to consider the effect of magnetically driven warping on the rotation of T 
Tauri stars. The projected rotation velocity of CTTS with masses M ^ 1 Mq is about 20 km s~^, 
only 10% of the breakup speed (e.g., Bertout 1989). This is at odds with the expectation that T 
Tauri stars are formed by the gravitational collapse of rotating molecular cores and the presence 
of disks surrounding the stars. Theories which explain the slow rotations generally invoke the 
interaction between the disk and the stellar magnetic field of a few kG (e.g., Konigl 1991; Cameron 
& Campbeh 1993; Shu et al. 1994; Yi 1995; Armitage & Clarke 1996). Since the growth time 
for disk warping is short, we may expect the inner disk of T Tauri stars to wander around the 
"preferred" perpendicular state. The star therefore experiences both spin-up and spin-down during 
its evolution, analogous to the behavior of X-ray pulsars (see §6.1). The net, secular spin-up rate 
is expected to be much smaller than that based on the canonical spin-up torque (~ M^JGMrm)- 



7. CONCLUSION 

In this paper we have identified a new magnetically driven warping instability which occurs 
in the inner accretion disk of a magnetized star (§2 and §4). Despite the uncertainties in our 
understanding of the magnetosphere-disk interactions (particularly the global magnetic field 
structure), the existence of the instability seems robust, and requires that some vertical field lines 
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(either from the star or intrinsic to the disk) thread the disk and get twisted by the disk rotation. 
The general consequence of the instabihty is that the normal vector of the inner disk (near the 
magnetosphere) will be tilted with respect to the stellar spin axis. We have also shown that 
the disk can be driven magnetically into precession around the spin axis due to the interaction 
between the screening surface current and the stellar magnetic field (§2). In addition, certain 
regions of the disk are subjected to resonant magnetic forces which may affect the structure and 
dynamics of the disk (§3). These magnetic effects on accretion disks have largely been overlooked 
in previous studies of accretion onto magnetic stars P^ . 

We have applied our theory to several different types of astrophysical systems, including X-ray 
pulsars, low-mass X-ray binaries, and T Tauri stars (§6). These applications should be considered 
preliminary, but they indicate that the magnetically driven disk warping and precession can 
potentially play an important role in determining the observational behaviors of these systems. Of 
particular interest is that the magnetically warped inner disk may provide a natural explanation 
for the longstanding puzzle of torque reversal as observed in a number of X-ray pulsars. Also, 
the tilted disk may be responsible for the rich phenomenology of time variability (such as QPOs) 
observed in weakly magnetized accreting neutron stars in LMXBs. 

Much work is needed to understand better the effects studied in this paper and their 
observational manifestations in different astrophysical systems. In our analysis, we have 
intentionally avoided (or bypassed), by using parametrized models, the uncertainties associated 
with magnetosphere - disk interactions (see, e.g.. Appendix A,B), but clearly the study of disk 
warping and precession based on more specific models (with or without outflows) would be 
useful. The role of intrinsic disk field needs to be examined further (see §2.4). There remains 
uncertainty in the description of nonlinear warped disks, and a full numerical simulation of the 
nonlinear development of the warping instability would be valuable. In the case of low-mass X-ray 
binaries, the effects of complex magnetic field topology (other than dipole) and general relativity 
should be included to access the QPO phenomenology (see §6.2). More detailed comparison with 
observational data will be useful. The role of magnetically driven resonances need to be studied 
further to determine whether they will produce any observable features. We hope to address some 
of these issues in the future. 

I thank Richard Lovelace, Phil Maloney, Dimitrios Psaltis, Marten van Kerkwijk and Ethan 
Vishniac for useful discussion/comment. I also thank Brad Hansen for informing me of some 
important references. I became puzzled by the spin-up/spin-down behavior of X-ray pulsars 
several years ago at Caltech through discussions with Lars Bildsten, Deepto Chakrabarty, Rob 
Nelson and Tom Prince. I acknowledge the support from a NASA ATP grant and a research 



^^We note that the effects studied in this paper are quite different from the situation considered by Agapitou et 
al. (1997), who studied the bending instability in the inner (sub-Keplerian) disk which corotates with the star (with 
the magnetic axis ahgned with the spin axis). Such a disk may or may not exist (see Spruit & Taam 1990). 
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A. MAGNETIZED ACCRETION DISK 

For an accretion disk threaded by a vertical magnetic field Bz, the steady-state angular 
momentum equation reads: 



Vr- = - 

dr S 



1 d f^3^^^^\ ^ BzAB^ 



dr \ dr J 2tt 



(Al) 



where I is the angular momentum per unit mass, S is the surface density, /S.B^ is the (/)-component 
of the magnetic field produced by twisting B^ and evaluated at the upper disk surface. Integrating 
over r and using the mass continuity equation M = —2TrrT,Vr, we obtain the conservation equation 
for angular momentum: 

Mlo = Ml + 2TTUir^m' + MIb, (A2) 
where Iq is a constant, 0' = di^/dr, and 



MIb = - drr'^BzAB^. (A3) 

Jr 

Equation ( [A2[ ) says that the rate of net angular momentum transfer through the disk, MIq, is 



equal to the sum of the rates of advective, viscous, and magnetic transport. We can rewrite ( A2) 
as 

(-) 

where 

J = (A5) 

The radial velocity is 

Vr = ^ (A6) 



For a Keplerian flow, n = y^GM/r^, I = ^GMr, and eqs. (0)-(|A|) reduce to (|^). The 
standard thin disk equations are recovered if we set Ib = 0. 

The above equations are quite general, but the actual expression for Ib depends on the 
behavior of the disk magnetic field, and thus should be viewed as being uncertain. Consider a 
specific ansatz: 

Bz = Bo(^) , AB^ = -CBz (A7) 



(see eqs. [ 2-27 l-[ p^^28| ), where Bq measures the magnetic field at the stellar surface (r = R). 



Assuming that a constant, we find 

" 3Mr3 3 



lB = ^ = lb%(^)\ (A8) 



- 29 - 



where Ir = \/GMR and = (^B^E^ / [MIr) is dimensionless. More general expressions (including 
spin dependence, multipole field and general relativistic effect) can be found in Lai (1998). The 
inner edge of the disk, r^, can be formally defined as where where d{l + lB)/dr = 0; we find 



R{2b' 



2 A 2/7 



\mVgmj 



(cf. eq. Il^l]). We then have 



^0 = K'^m) + h{rm.) 







and 



J 



6 V r 



1/2 



7-1 ^ 



(A9) 



(AlO) 



(All) 



This would suggest ^7 = at r = r^- But of course as r approaches the inner disk boundary, this 
expression breaks down, as the inertial term of the radial equation must be taken into account 
(see Lai 1998 for a magnetic slim disk model). 



B. MAGNETOSPHERE RADIUS FOR GENERAL FIELD GEOMETRY 

In the main text, we have intentionally avoided the precise definition of the magnetosphere 
radius for general orientations of /, u) and (i. This is because the exact determination of 
requires a detailed model of the structure of magnetosphere-disk boundary layer, for which no 
definitive theory exists (useful discussions are contained in the references cited in §1). We offer the 
following three possibilities. They should be viewed only as an educated guess, although they all 
give a scaling relation as in eq. (1-1). 



(i) Ansatz 1: If the static vertical field B^ = — {fi / r^) cos f3 cos 6 threads the disk (see §2.3), 
it will affect the angular momentum transport in the disk through the magnetic stress oc B^AB^j,, 
where B^ = —C,Bz is the field created by winding Bz (see eqs. [1A^ - |1A3|| ). When the condition 



+ r^BzAB^ = (Bl) 

dr 

(where I = VGM r is the angular momentum per unit mass) is satisfied, no viscosity is needed to 
induce accretion. We may thus use this condition to determine the inner edge of the Keplerian 
disk, giving 

1 /7 

r„=(2Ccos^/3c„s^«)'"('^) . (B2) 



(ii) Ansatz 2: Here we shall follow the consideration similar to that in Arons (1993). Assume 
that the disk is largely diamagnetic, but strong dissipation exists at the boundary layer, where 
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plasma stress and magnetic stress balance, i.e., 



pVrVs 



4tt 



The vertical magnetic field at the disk inner edge {r = rm) is given by 



Bz{rm) 



4/i / r. 



1/2 



cosx 



(B3) 



(B4) 



(Aly 1980). Shear in the boundary layer induces Ai3^ = —CbBz with ^ 1- Using = GMj r 
and Vr = -M /{A-KvpH), we find 



8Cb 2 
— cos X 



2/7 



1/7 



GMW^ 



(B5) 



Since cosx = cos/? cos — sin/3sin0sina;t, eq. ( |B5| ) implies that the magnetosphere boundary is 
modulated at the spin frequency w. If we average the magnetic stress over time, we may replace 
cos^X by [cos^/3 005^6* + (1/2) sin^/Jsin^^]. 

(iii) Ansatz 3: This is similar to (ii), except that we now assume that the static vertical field 
threads the disk [see eq. ( ^-25 )], so that at r = r^, 



Bz{rm) = — T' "^os 6 H ^ — ^ I sin /3 sin 9 sinujt. 



(B6) 



The condition ( p3|) then gives 



8Cb\ 
7r2 j 



2/7 



sin P sin sin ivt 



vr /2if\i/2 
— ( 



4 Vr„ 



cos /3 cos 9 



14/7/ 4 \ 1/7 
( ^ . ) 



(B7) 



We see that in general, the inclined rotating dipole will produce a time- varying magnetosphere 
boundary. This variability may be in addition to that associated the winding-up and reconnection 
of field lines (see the end of §2.4). 



C. MOTION OF DIAMAGNETIC BLOBS 

It has recently been suggested that under certain conditions, the accretion disk onto 
magnetized object may consist of diamagnetic blobs (King 1993; Vietri &; Stella 1998). Here we 
study how the stellar magnetic field affects the motion of diamagnetic blobs which lie in circular 
orbits of the disk plane. Each blob moves independent of each other, and is subjected to the 
gravitational force from the central star and the magnetic drag force which arises when it moves 
across the field lines (see below). It is not clear that a realistic disk will behave as a collection of 
individual blobs, nor is it clear that the blob can survive for a long time (e.g., the blobs may be 
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subjected further Kelvin-Helmhotz instability which tends to break up the blobs) — These issues 
are not addressed in this paper. Our purpose here is to understand the dynamics of the blob 
under the afore-mentioned assumptions. 

Vietri and Stella (1998) has studied some aspects of this problem. They assumed that the spin 
axis is aligned with the angular momentum axis of the orbiting blob. They showed that vertical 
resonances exist in the region where iojl < < 3uj/2 (where uj, fl are the spin frequency and 
orbital frequency), and suggested that the resonances can pump the blobs out of the equatorial 
plane. Here we consider the more general case where the orbital angular momentum of the blob is 
misaligned with the spin axis. We show that even without the resonances, there exists a magnetic 
torque which can induce tilt of the orbit of the blob. 

Our setup is given in Fig. 1. The magnetic drag force on the blob results from the Lorentz 
force on the screening current on the blob's surface. The blob loses energy by exciting Alfven 
waves in the surrounding plasma. The drag force per unit mass is given by 



and the drag time scale is 



drag 



AVj 

Td 



CAm 



(CI) 



(C2) 



Here ca is the Alfven speed in the interblob plasma (where the field strength is B), and m and d 

is the mass and characteristic size of the blob (Drell, Foley & Ruderman 1965). We shall assume 
that the interblob field is that of the stellar dipole (thus neglecting possible screening due to the 
blobs): 



B. 



^sinxcos(0- ^^), 
^sinxsin(^ - ^^), 

3 COSX. 



The relative velocity between the blob and the field line [at location (r, (t),z = 0)] is 



AV = Qrcl) — Qs<^ X r. 



(C3) 

(C4) 
(C5) 

(C6) 



The projected relative velocity perpendicular to the field line is then 

(B • AV) B 



AVi 



AV- 



|B|' 



AQr (j) + Ljr sin (3 cos (j) z 
T r 

An sin X sin(^ — ^^) — u cos x sin (3 cos (j) 



C 



2 sin X cos(^ — (j)^) f + sin x sin(^ 



cosx^ 



(C7) 
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where 

C = l + Ssin^xcos^l^-^M)' (C8) 

Ai7 = O-wcos/?. (C9) 

The drag constant can be written as 



(4tw)'", (Cll) 



where 

(po is the density of the interblob medium). The second equahty of eq. ( piO| ) is vahd only when 
CA is less than the speed of light. 



C.l. Magnetic Torques 

The magnetic torque on the blob (per unit mass) is given by 
N = 



2 r 

(An z - uj sin P cos (j) <j) 



r 

TdO 



AQ sinxsin( 



to cos X sin P cos ( 



sm X sin[ 



z + cos X ' 



(C12) 



To average over 0, we shall approximate C^/^ by and C~^/^ by C2 , where C\ and C2 are 
constants in the range of 1 — 4. Since the Taylor expansion of C^^^ or C~^/^ contains cos 2n( 
(where n is an integer), we can be sure that no new term (with different dependence on the 
would appear if the exact expression of or 

were adopted. We find 



angles) 



C^^ [ An Z--UJ sin py 



1 



-a 



-1/2 



TdO ^ \ 2 J 2TdO 

+ (^An sin X cos x sin + uj sin /? cos^x) V 
— ( An siv?x + ^ sin /3 sin x cos x sin ( 



An sin X cos x cos x 



Averaging over the rotation period and using the identities ( p-9| )-(2-ll) we obtain 

„2 



((N)) 



-— C^^ (Anz--ujsinPy 
TdO \ 2 



(C13) 
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-C. 



-1/2, 



w cos^O + - (Q cos (3 — Lo) [?> cos^O — 1 



sin f3 y 



An sin^e + -n sm^/? (3 cos26' - 1 



(C14) 



In the case of /? = (i.e., aligned ui and /), the only nonzero torque is along the z-axis: 
.2 / 1 



((N)) 



a 



1/2 



(C15) 



This is simply the magnetic breaking torque which tends to drive the blob toward corotation with 
the magnetic field lines. For /? 7^ 0, there is another torque in the y-direction, which tends to tilt 
the orbit of the blob. In the absence of other forces, the tilt angle /? evloves according to 



1 



dp_ 

dt 4Tdo 



sin P 



a 



'^/^ (3cos26'-l)cos/3 + 



/2 



sin^6l - 2C 



1/2 
1 



(C16) 



When a; <C ri, this equation has a similar structure as eq. ( ^-22 ). Clearly, under certain conditions 
(when the quantity inside the sqare bracket is positive), there is an instability where P tends to 
grow toward the perpendicular state (/? = 90°). The growth time is of order t^q. 



C.2. Resonances 

The vertical magnetic force (per unit mass) on the blob is given by 
Fz = —T^^C^^'^Lor sin (3 cos (j) 



-^do^ ^/^rcosx 
The equation of motion for the vertical motion is simply 



AOsinxsin((/) — 0^) — a;cosxsin/3cos(/) . (C17) 



+ niz = Fz. (CI8) 



First consider the case where /3 = 0. The time-dependence of the force is as in 
C~i/^ sin((/) — ujt), which can be written as a sum of sin(2n + 1)(0 — ojt). Thus the resonance 
conditions are (2n + l)(u; — 17) = ±17^, or, for nz = n (Keplerian disk): 

2n 2n + 2 

""^2^^' °' ""^2^^' (n = 0,l,2,---). (C19) 

Thus the resonance "band" lies inLj/2 < n < 3uj/2 (but excluding uj = n). This was first 
identified by Vietri & Stella (1998). Note that for uj = 0, the resonance = $7 is always satisfied, 
which implies that the blob will be driven out of the equatorial plane. 
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Now for the general (3^0 cases. If we treat C as being independent of time, then it is easy 
to identify the following resonances: (i) fi^ = J7 (which is always satisfied for a Keplerian disk); 
(ii) u!±Q, = ifi^; (iii) 2u!±Q = ifi^. To take account of the time-dependence of C, we note that 



C= 1 + 3 



(3 (3 
sin 9 s\v? — cos + uit) + sin cos^ — cos ( 



ujt) + cos 9 sin (3 sin (j) 



(C20) 



Thus C^/^ or C^l'^ can be written as a sum of cos(2n^ it mwt) and sin(277,^ it mwt) (where 
n, m = 0, 1,2, • • •). We can then show that resonances occur when 2nS7 = mu (for Jl^ = f2), 
i.e., all blobs (with different orbital radii) are subjected to resonant forcing. We expect that the 
high-order resonances (those with large n and m) are weak, although we have not tried to quantify 
them. The ubiquity of the vertical resonances implies that the disk consisting of diamagnetic blobs 
tends to thicken due to the magnetic drag force. 



The radial force (per unit mass) on the blob is given by 

Fj. = 2T^C~^^'^r sin x cos{(j) — (f)^) Ail sin % sin(^ — (j)^) — lo cosx sin/3 cos ( 



(C21) 



Similar consideration reveals the existence of a large number of epicyclic resonances in the disk. 
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Fig. 1. — Coordinate system used in the calculation. The z-axis is along I, the angular momentum 
of the disk (or a ring of disk). The angular velocity vector of the star, u) = uu (where u is the 
unit vector) is inclined at an angle /? with respect to I, and lies in the yz-planc. The stellar dipole 
moment fi = fifi (where fi is the unit vector) rotates around uj, and the angle of obliquity is 9. In 
the cartesian coodinate, we have fi = sin x cos (/)^ x + sin x sin (j)^y + cos x -2- The axis z' (along the 
spin axis) is used for studying disk warping, and u± is the unit vector perpendicular to uj and lies 
in the yz-plane. 



